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Abstract. We present results of an analytic and numerical calculation that studies 
the relationship between the time of initial foot contact and the ground reaction force 
of human gait and explores the dynamic principle of center of mass. Assuming the 
ground reaction force of both feet to be the same in the same phase of a stride cycle, 
we establish the relationships between the time of initial foot contact and the ground 
reaction force, acceleration, velocity, displacement and average kinetic energy of center 
of mass. We employ the dispersion to analyze the effect of the time of the initial foot 
contact that imposes upon these physical quantities. Our study reveals that when the 
time of one foot's initial contact falls right in the middle of the other foot's stride cycle, 
these physical quantities reach extrema. An action function has been identified as the 
dispersion of the physical quantities and optimized analysis used to prove the least- 
action principle in gait. In addition to being very significant to the research domains 
such as clinical diagnosis, biped robot's gait control, the exploration of this principle 
can simplify our understanding of the basic properties of gait. 



PACS numbers: 87.85. G, 87.85.gj, 87.55.de. 
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1. Introduction 

Gait analysis plays an important role in exploring laws of human motion by gait 
parameters via biomechanical methods. Many studies have shown that gait parameters 
are significantly symmetric [11 [2] and can be understood in terms of segmental 
kinematical and kinetic physical quantities, while walking [SI H]. It has been shown, 
[3 El [TJ E] that by developing gait parameters into evaluation indexes, one can assess 
the causal relationship between physical injury and gait ability, which in turn has been 
applied in rehabilitation therapy with success [9l El [TOl [TTl [121 1T3] . 

The essence of human center of mass (COM) motion is actually the periodical 
change which is under the effect of external forces such as, ground reaction force (GRF), 
gravity and resistance [H]. Previous studies indicate that the law of COM motion is a 
more reliable gait evaluation method [I5]. The common features of the frequently- used 
gait indexes lie in the fact that normal human gait parameters have been regarded as 
the criteria to evaluate rehabilitation. However, the incomplete symmetry of human 
shape, coordinate and strength has formed the uniqueness of human gait [17]. This has 
brought difficulty to the establishment of standard gait parameter index. Consequently, 
exploring a principle that is relevant to gait parameters such as cycle time, cadence 
or stride length in a normal human gait has become essential to the study of gait 
biomechanics [H]. There are a number of issues we wish to explore in our calculation. 
An important question is whether the relationships between the GRF, velocity of COM, 
average kinetic energy and the time of initial foot contact (tIFC) in a stride cycle can be 
established. In particular we wish to examine if and to what extent such relationships 
explore the principle behind the gait characteristics. In the present work, we wish to 
advance an attempt towards an approach where the precision of the calculation reaches 
new levels. 

The discovery of many basic principles originated from the study of animal 
movements [H]. By examining the GRF acted on the normal human gait (bare- footed) , 
we establish the relationships between the GRF, velocity of COM, average kinetic energy 
and the tIFC in a stride cycle. In this contribution we present our recent results on 
kinetic regularities of COM. Another crucial issue in this attempt concerns the prediction 
of action in gait. Using analytic and numerical techniques, we shall identify the action 
in gait and explore the possibility of the least-action principle in gait. The rest of the 
paper is organized as follows: 

In Sec. II we describe our procedure to establish the relationships between the 
GRF and the tIFC and discuss the aspects of determination of the force distribution 
in anteroposterior, mediolateral and longitudinal directions. We address the problem 
of the force distribution in transverse, sagittal and frontal planes and optimization of 
dispersion of force. We present a description of the determination of working parameters, 
such as acceleration, velocity and displacement of COM and establish the relationships 
between the kinetic regularities and the tIFC in this section. We conclude this section 
by addressing the effect of the tIFC upon the average kinetic energy of COM. Our main 
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results are presented and discussed in Sec. III. Finally, we present our conclusions in 
Sec. IV. 

2. Analysis Strategy - Forward Dynamic Method 

2.1. The Effect of tIFC upon GRF 

The analysis of COM's dynamic characteristics to evaluate the gait features is rather 
non-trivial. The inverse dynamic method [IHl UHl [16] has given mixed results with 
systematic and statistical errors as major sources of uncertainties in the analysis of COM 
dynamics characteristics [20| [2Tj. An attempt to indicate the dynamics characteristics 
of the whole body COM by one particular segment of body will highly underestimate 
the results. The target of this work is a calculation that gives more precise estimates 
of dynamical measurements in question. We use forward dynamics method [221 EH] to 
illustrate on how tIFC determines the force distribution and establish a relationship 
between the tIFC and its force dispersion. 

Since dynamic characteristics of COM in gait is governed by kinematics and kinetics 
of COM, therefore, for constant gravity, the analysis of force upon COM will compliment 
the analysis of GRF. Also, GRF is caused by the body segmental movement driven by 
the transarticular muscles and eventually by the contact of the foot to the ground, so 
one foot's GRF includes the longitudinal GRF, frontal friction and sagittal friction, i.e., 

F(t) = F,(t)i + F,(t)j + F,(t)k, (1) 

where F{t) is the GRF of one moment at a stride cycle and Fx[t), Fy{t) and Fz{t) 
represent the components in three directions, respectively. If Flit) and F^it) {i = x, y, z) 
represent both feet's sagittal and frontal frictions and longitudinal GRF, respectively, 
then the above equation shows a similarity of both feet's GRF distribution when 
expressed in biped gait features [T8j. Consequently, we assume that in each stride cycle, 
the left and right foot have an identical distribution of GRF in all three directions. This 
would imply that the GRF variations are only caused by the tIFC of both feet. Setting 
T as one foot's stride cycle time, the initial phase of one foot is equal to zero and that 
of the other foot to (i.e. tIFC), then corresponding to one foot's GRFs (for example left 
foot) 

Fi{t)=F,{t), Fl{t) = Fy{t), F[{t) = FM, 
those of the other foot are 

F;(t) = F,.(t + t„), Fl{t) = Fy{t + to), F:{t) = F,{t + to) 
and Eq. ([1]) becomes 

F(t, to) = F^{t, to)i + Fy{t, to)} + F,(t, to)k. (2) 

Since gait is a continuous and periodic movement, therefore while walking at steady 
speeds, -F(t, to) is the GRF when tIFC is to and the stride cycle time is t. If 
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F{t,to) = Fit + nT.to + uT) (n = 1,2,3. . .) holds, then the GRF and gravity {W) 
in a stride cycle that act on the impulse of COM will follow the following equation: 

/ = 4i + lyi + /,k = 0. (3) 

To calculate the impulse characteristics of each foot's GRF in gait when both feet have 
same GRF, we analyze two phases of one foot's stride cycle. Let and represent the 
stance phase and swing phase, respectively, then in the anterposterior direction, Eqs. 
(E]) and ^ yield 

f F^{t)dt = and [ F^{t)dt = f F^{t + to)dt. 
Jo Jo Jo 

Since T = + and j^"^'^'" Fx{t)dt = 0, therefore 



pis ris 

/ F,{t)dt= / F^{t + to)dt = 0, 
Jo Jo 



whereas in the mediplateral direction, we get 

I Fy{t)dt = / Fy{t + to)dt. 
Jo Jo 

On the other hand, the calculated contribution in the longitudinal direction gives 

F,{t)dt = [ F,{t + to)dt = ]- I Wdt. 
Jo ^ Jo 

The above equations give the impulse characteristics of each foot's GRF in gait. 



2.2. Dispersion of GRF 

To analyze the dynamics that emerge in the GRF distribution at to = 0,T/2 and T, 
we bring forward the concept of dispersion of GRF. If cTj (i = x, y, z) and Fi denote the 
GRF dispersions and average values of GRF in three direction, then the corresponding 
correlation between the force dispersion and average GRF can be written as 



T T 
t=0 t=0 



{Fi{t,to)-mo)f/Tj2^^ (4) 



where Fi{t,tQ) is the same as that in Eq . ([2]) and Ylt=o ^ number of At within 

the range [0,T]. For sequential values of to within the range [0,T], we calculate the 
values of ai and use them to evaluate the effect of tlFC upon GRF, velocity, position 
and kinetic energy of COM. We use the average value of GRF from 20 subjects to check 
the reliability and accuracy of this method. This signature is confirmed in Sec. Ill of 
this study. 
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2.3. Regularities of velocity and displacement of COM 

Ground reaction forc^ is the result of body segments action on the ground via foot 
and determines the kinetic regularities such as acceleration, velocity and displacement 
of human body COM. Knowing the GRF and weight, the acceleration of COM at a 
given instant in a stride cycle using Eq. (jj]) can be written as 



Expressing the acceleration of COM in component form, Eqs. ([T]) and ([5]) reveal that 
the tIFC and GRF have the same effect on the acceleration of COM. 

The absolute motion of human COM relative to absolute inertia reference frame 
is composed of convected motion and relative motion. Following Kokshenev [11], we 
define gait in accordance with motion in Eq. ([3]), as "walking at steady speeds", thus 
we regard the convected motion a constant parameter. Using Eq. ([5]) at tIFC = to, the 
COM velocity in a stride cycle is given by 



where fo(to) is the initial velocity of COM in relative motion at the beginning of a stride 
cycle. It seems that we cannot confirm the magnitude of v{t,to) in Eq. ([6]) by kinetic 
method (since everyone's gait speed is different). However, the distinct feature is that 
once GRF and tIFC are determined and when it observes the motion in Eq. ([3]), Vo(to) 
must have a unique solution. We can then establish the relationship between the initial 
velocity and acceleration in relative motion in a stride cycle 



Eq. reveals that v{T,to) = v(to), which implies that the cycle of velocity of COM 
is in accordance of v{T + nT,to + nT){n = 1,2, • ■ •), that is to say, the end of one 
stride cycle marks the beginning of the next stride cycle. This explains the periodical 
characteristics of velocity of COM in gait. What needs to be further illustrated is that 
the initial velocity calculated using Eq. ([7]) refers that at the beginning of a stride cycle, 
the body is in the state of steady speeds. From Eqs. ([5]) and ([7j) we conclude that 
tIFC has determined the initial velocity in relative motion, which has nothing to do the 
convected velocity. 

Just like the velocity of COM, displacement of COM also involves absolute, 
convected and relative displacements. We define the COM initial displacement in 
relative motion of different tIFC as So(to)- Using Eq. ([6]), the displacement of COM in 
relative motion at any moment has the following form: 



a{t,to) = F{t,to) -Wk. 



(5) 




(6) 




(7) 




(8) 



I Wc assume a constant gravity and a negligible air resistance. 
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and the correlation between the initial displacement of COM and velocity of COM in 
relative motion in a stride cycle can be written as 



Similar to the behaviour of the initial velocity in relative motion, we verified that the 
initial displacement seems to be independent of gait velocity, cadence and stride length. 

2.4- Average kinetic energy 

Nature has always minimized certain important quantities when a physical process takes 
place [21]. Bipedal walking has enabled the continuous evolution of human gait [25] and 
eventually it has brought about the optimized gait [26] and thus became a behavioral 
trait. In this behavior, the force, acceleration and velocity all have their minimal 
values when to = \T. To understand the physical significance of these gait dynamic 
characteristics we explore the issue of energy consumption. In gait, human segment 
movement is a combination of agonist, antagonist and synergist. The movements such 
as stretch or flexure all consume mechanical energy. On the other hand, E'fc > in COM 
kinetic energy, ^^-E'fc(^) > in total kinetic energy and Ek{t) > in average kinetic 
energy in a stride cycle, whereas the corresponding potential energy counterparts are 
zero respectively. Using the COM kinetic energy in relative motion, the description of 
mechanical energy consumption in gait simplifies to 



and can be easily resolved in component form. Knowing iJj, we can set up the 
relationships between tIFC and COM kinetic energy by defining t and to in the interval 



2.5. Experimental Details 

The experimental measurements were carried on the following set of equipments: Simi 
Motion 7.0 Three-Dimensional Movement Analysis System; three Kistler 40 x 60 [cm?) 
force plates; force plate frequency: 2000iJ2;, with a systematic uncertainty of ±1%. The 
assembled force plate position has been fixed by gradienter to ensure the force plates are 
on the same plane. The measurements are taken at the sample frequency of lOOOif^;. 
Before each measurement, the equipment is examined and returned to zero-position. 

Twenty female subjects with mean age 20.63±0.76 years, mean height 163.12±3.72 
cm and mean weight 45.74 ±3.30 kg participated in the study. A few trials of each level 
gait item were administered to subjects as they ambulated on instrumented positions. 
All the patients agreed to participate in the research, and signed freely an informed 
consent form and study was carried out according to the existing rules and regulations 
of our institute's Ethnic Committee. Before the test, all the subjects were thoroughly 
briefed about the procedures and matters needing attention so that they all understand 




(9) 




(10) 



[o,r]. 
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the purpose and the requirements of the test. Each subject's medical history is inquired 
so as to exclude subjects with diseases such as pathological change, deformity or injury 
to make sure that their physical conditions would meet the requirements of the test. 
When measuring their gaits, we start from the subjects' standing position, bare-footed 
(both feet disinfected by 75% of ethanol). A pre-test is to guarantee that after they 
walked three steps, the subjects step on the platform and to make sure that each force 
plate can measure one foot's stance phase data separately. When the subject's gait is 
found to be quite abnormal, for example, it is obviously discontinuous, she would be 
asked to perform again so that the recorded data meet the requirements of the test. 

Based on the fact that the longitudinal GRF in gait is apparently greater than the 
GRFs of the anteroposterior or mediolateral direction, we take the signal of longitudinal 
GRF to identify the instants of initial foot contact and of terminal stance. The gait we 
study is walking at steady speeds; therefore, we de-noise the signals of the platform and 
examine the test signals would meet the requirements in the predetermined accuracy 
range. We rate the gait cycle time by percentage, normalize the weight and standardize 
the GRF from three directions. We count on average of the processed data of 20 subjects' 
three-direction GRF from individual stride cycle to conduct our research. 

3. Results and Discussion 

3.1. Spatial GRF, COM velocity and displacement 

To explore the spatial GRF, we examine Eq. ([2]) by using the GRF numeric of 20 
subjects. Fig. [1] collects and displays the results of anterposterior, mediolateral and 
longitudinal GRF. The GRF has been standardized {Fi(t)/W)) and stride cycle time 
is rated by percentage. Fig. [HA. - C shows that tIFC has an effect upon the GRF in 
three directions. Using three sets of component data from GRF of the subject's gait, we 
obtained curved surface effective plots and the quantitative relationship between tIFC 
and GRF resultant forces on three planes. Fig. [Tp - F indicates that no matter what 
changes to undergoes, GRF resultant forces in these three planes remain to be closed 
curves. This significantly compliments and confirms the dynamics behind Eq. ([3]), i.e., 
when the distribution of one foot's GRF is determined and both feet's GRFs remain 
the same, the momentum of COM always remains unaffected. 

Fig. [Tp - F also reveals the geometrical characteristics of the plane resultant force. 
On the transverse plane, when t^, = 0, T, the gait becomes a jump and its frontal 
resultant force becomes a straight line. The resultant force changes into a symmetrical 
butterfly about F^, which could imply a possible transformation at t^ = \T. On the 
sagittal plane, when t^ = 0,T, the resultant force Fr^z forms the longest closed curve 
(the length of which is calculated by § Fxz{t,to)dFxdFy) and more or less shortest at 
to = \T. On the frontal plane, when the resultant force develops into a straight line at 
to = 0,T, the resultant force becomes a symmetrical butterfly about F^ at to = |T. 

Having developed the relationship between tIFC and GRF's resultant force on three 
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Figure 1. Schematic illustration of the observed spectrum for tIFC and GRF showing 
the relationship between tIFC and GRF in three directions (Fig.lA - C), together with 
tIFC and GRF resultant force in three planes (Fig. ID - F), compared with a typical 
tIFC and dispersion of force correspondence (Fig.lG - I). 

planes, we calculate the GRF dispersion in three directions using Eq. To confirm 
the results beyond doubt, we identify the GRF dispersion by analyzing the effective 
plots shown in Fig. [in - I. It is clear from the plots that at to = ^T, ax has its 
maximum value at the bottom region of "W" shape, ay has a global maxima and o"^ 
a global minima. Using the numerical values, min^ax) = 0.0695, max{ax) = 0.1319, 
min{ay) = 0.0000, max(aj,) = 0.0355, mm(or^) = 0.1086 andmax{a,) = 0.9144, of GRF 
dispersion obtained in this case study, we find that the largest differences between the 
maximal and minimal values are 0.0624, 0.0355 and 0.8058, respectively. This means 
that the dominant effect to the global GRF is the longitudinal GRF. The minimal 
GRF dispersion on longitudinal direction occurs at to = |T. The analysis of the 
distribution of GRF in three directions indicates that when to = |r, the anteroposterior 
and longitudinal dispersions of GRF are almost the least and the mediolateral one is the 
largest. The closed curve shaped by the two components on the sagittal plane seems to 
be the shortest. The GRF resultant force on the transverse and sagittal plane is shown 
as symmetric butterfly. 

Turning our attention towards the kinetic regularities, we plot the effect of tIFC 
on COM velocity in Fig. O We can see two platforms in Fig. |2]A. - C, which emerge 
at the two extremes of tIFC rated by percentage forming a concave region around tIFC 
whereas the velocity describes a convex region around tIFC. The COM velocity in the 
longitudinal direction indicates that tIFC entails the vertical velocity of COM to have 
a distribution of tIFC in the form of a saddle. In order to further understand the effect 
of tIFC exerting on the velocity of COM, we analyze the variations of the velocity of 
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Figure 2. Effective plots showing relationship between tIFC and COM velocity in 
three directions (Fig.2A - C), tIFC and COM velocity in three planes (Fig. 2D - F) and 
tIFC and dispersion of COM velocity (Fig.2G - I). 

the direction of COM velocity is a straight line at to = and T and the plane velocity 
becomes a symmetrical closed curve about v^^ axis at to = \T. In the sagittal plane 
the closed curve formed by the plane velocity is the longest (the length of the closed 
curve is calculated by § v{t,to)dVxdVz) at to = 0,T, and is shorter (not the shortest) 
at to = and T. In the frontal plane, when to = 0,T, the velocity becomes a straight 
line and the resultant velocity becomes a symmetrical closed curve about axis at 
to = \T. Having obtained the estimated value of COM velocity, we are ready to set 
up the dispersion of relationship between tIFC and COM velocity (as shown in Fig. 
[2n - I) using Eq. We notice that at to = |T, a^- and cxy have the global minima 
and has the global maxima. The estimates, min{ax) = 0.0541, max{ax) = 0.1652, 
min{ay) = 0.0000, max{o-y) = 0.0510, min[az) = 0.0549 and max{ax) = 1.3502 in our 
case study, result in largest differences of maximal and minimal values of 0.1111, 0.0510, 
and 1.2953 respectively. This means that the dominant effect upon COM velocity lies in 
the longitudinal and anteroposterior directions, where the dispersions have the minimal 
values at to = \T. 

The kinematic regularity of the displacement of COM in three directions is shown 
in Fig. I3JA. - C. We see a geometric distribution of convex corners indicating that tIFC 
changes the position of COM in anteroposterior and mediolateral direction. The bimodal 
plot indicates the complexity of longitudinal direction exerted by tIFC to the position of 
COM. The estimated values of the COM displacement components are used to evaluate 
the values and directions of COM displacement in three planes. To further analyze the 
effect of tIFC upon the position of COM, we explore the changes of COM positions from 
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Figure 3. COM displacement as a function of tIFC. The directional and planer 
components arc shown in Fig.SA - C and Fig. 3D - F, respectively. Fig.SG - I shows 
the dispersion of COM displacement. 

These effective plots show that in the transverse plane the COM displacement is 
a straight line at to = 0, T. On the other hand a closed curve is formed by and Sy 
and is symmetric about Sx axis at to = |T. On the sagittal plane, when to = 0,T, the 
closed curve shaped by the sagittal displacement is the longest (the length of the closed 
curve is calculated by j s(t,to)dSxdyz); when to = |T, the curve is rather shorter (but 
not the shortest). On the frontal plane the frontal displacement becomes a straight line 
at to = 0, T and a symmetric closed curve about Sy axis 

Using the numerical values, min{ax) = 0.0042, max{ax) = 0.0252, min{ay) = 
0.0000, max{ay) = 0.0081, min{az) = 0.0039 and max((j^) = 0.2147, obtained in this 
case study, we find that of the GRF dispersion ax and Cy have a global minima and o"^ 
has a global maxima at around to = ^T, again resulting in largest differences between the 
minimal and maximal values with a magnitude of the order 0.0210, 0.0081 and 0.2108, 
respectively. This implies that the greatest effect to COM displacement comes from 
components on longitudinal and anteroposterior directions while the minimal dispersion 
of displacement on the longitudinal and anteroposterior direction emerges at to = |T 
(See Fig. [3p - I). 

Accordingly, at to = ^T, COM velocity in relative motion presents its symmetric 
closed curve in transverse and frontal plane and the curve length in the sagittal plane is 
approximately the minimal. In a stride cycle, the dispersion of COM velocity in relative 
motion has the global minima on the anteroposterior and longitudinal direction while 
the global maximal value on the mediolateral direction. The effect of tIFC upon the 
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COM velocity and displacement does not depend on gait velocity, cadence or stride 
length whereas the effect of the tIFC upon COM acceleration is in accordance with 
the effect it exerts upon the GRF. The COM acceleration and GRF show more or less 
identical behaviour. Fig. HJA. - C demonstrate the effect of tIFC upon three components 
of GRF, COM velocity and COM displacement on three direction. As is clear from Fig. 
HP - F that at half stride cycle {t = T/2), the dispersion of GRF and COM kinematics 
are minimum. In comparison to GRF and COM velocity, the COM displacement shows 
a sharp dip with least minimum. 




stride ivtle Stride cycle 



Figure 4. Relationship between tIFC and dispersion of force. 



3.2. Anteroposterior, Mediolateral and longitudinal average kinetic energies 

Finally, we show the COM kinetic energies as a function of stride cycle in anteroposterior, 
mediolateral and longitudinal channels in Fig. 0\ - C, together with corresponding 
transverse, sagittal and frontal energies in three planes (Fig. |5p - F). Fig. [5]G 
- I indicates that tIFC makes COM average kinetic energy the global minima on 
anteroposterior and mediolateral direction while that on the mediolateral direction a 
global maxima at to = ^T. Comparing the relationship between the tIFC's velocity 
of COM and the position of COM, tIFC contributes a symmetric distribution of 
kinetic energy of COM. Since min{Ex{t,to)) > max{Ey(t,to)) and max{Ey{t,to)) < 
min{Ez{t,to)), the anteroposterior and longitudinal average kinetic energies determine 
the COM average kinetic energy, when tlFC falls at the 20%, 50%, 60% and 100% of a 
stride cycle, the kinetic energy of COM in three planes varies rather considerably. On the 
transverse and frontal planes, the effect is similar whereas the sagittal plane undergoes 
a non-trivial change. We apply the dispersion of kinetic energy of COM to evaluate 
this effect. We notice, from Fig. |5]J, that COM average kinetic energy is minimum at 
to = |T, which implies a minima for the COM mechanical energy consumption and 
signatures the physical significance of gait. 
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Figure 5. Relationship between tIFC and kinetic energy. (Fig.5A - C tIFC and COM 
kinetic energy in three directions; Fig.5D - F -tIFC and COM kinetic energy in three 
planes; Fig.SG - I - tIFC and average kinetic energy. ) 



3.3. Least- action principle in gait 

The intriguing question that arises from the above discussion if there is a symmetry 
or a physical principle that is hidden in the phenomenon that makes the COM kinetic 
energy to reach its extremum at to = ^T. In order to explore this issue, we pursue the 
technique of optimizing the dispersion objective function |18] 

T T T 

+ /(^ + to) - HQ? = /w' + E /(^ + to) 



\2 



T T T 

+2 ^ f{t)f{t + 1,) - 2 ^(/(t) + f{t + to)) + a {to))' 5^ 1 (11) 
11 1 

The action of GRF and gravity are the reasons for human COM motion changes. Let's 
first analyze the GRF. The average value of the resultant force in the anteroposterior 
and mediolateral direction F^ito) and Fy{to) is zero and in the longitudinal direction 
Fz{to) is one (when the weight has been normalized). Therefore, for a known gait, 
Y.l{F^m. EiiFyit))^ Y.lmt))\ E^AF.it + to))\ EliFyit + to))' and EiiF.{t + 
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to))^ are all constants and Eq. fllip simplifies the dispersion optimization to 



T 



1 

T 



1 

T 

min^F^(to) = + (12) 

1 

To find a solution of the optimization problem, we need to set up a function of GRF 
that changes with time in three directions. It seems that the pulse periodic GRF fits 
the criteria and hence we use segment trigonometric functions 

F,{t) =asm{27i^), Fy{t) = bsm{n^) F,{t) = csm{7i^), 

-'■ s -'■ s -'■ s 

where a, b, c stand for the GRF after normalization in three directions and Tg is the 
stance time. For a small enough test frequency, the above equation transforms to 

Z"^" t t + to 

mmipF^ito) = A sin(27r— ) sin(27r ° )d^ 

Jo J^s 
/■^s t ^ -I- 

maxipFyito) = -B sin(7r— ) sin(7r ° )dt 

Jo -'■s -'s 

Z"^'' t t + to 

mm il: F Ato) =C sin(7r— ) sin(7r— -^)cit. (13) 

Jo -'s 



Let's take longitudinal GRF in Eq. (1131) as an example. In order to get the antiderivative 
of the integrand of integral variable t, we detach t and to of trigonometric function to 
obtain 

/ sin(;^7r)sin(^^^7r)c/t = / (sin2(^) cos(^) 

JO -'■s Jq Ig Ig 

+ sin(^)cos(^)sin(^))rft, (14) 

Now we can transform Eq. ( IT4|) into the integral of integral variable t while regarding 
the GRF as segment function. The transformation reduces the pulse periodic GRF 
functions into three piecewise functions at three intervali§|: [0,T^], [Tuj,Tg] and [Ts,T] 
which, to the confirmed GRF, a, b, c are constants and A, B, C (related to a, b, c) are 

§ While walking + T„ = T and T, > > 0. 
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also constants, yields the following form of longitudinal GRF: 



- smi TT— I - 
2 ^ T/ 

< to < 

1 . / to , 

- smlvr— I - 

2 ^ T/ 

^1 . .to- 
+ - sml — — 
2 ^ T, 



f — T f 

— TTl COSlTT— 

2T. ' ^ T, 



2T, 



-vr cos vr- 



T 



-TT - 



2T, 



T t 

-TT COS 



T 



(15) 



1 



< to < 



f —T 



sm 



2T< 



T t 

-U) N / ''O 

TT COS — 



T 



-vr 



The contribution f^^^^vr) sin(7r|¥-) = in the interval [0,T^„], thus giving the to 



which in turn yield (from Eq. 
contribution (^^^^vr) sin(7r|f ) + 



f fTSj) ) a maxima 



C7r^ 
2 > 



On the other hand, 



= 0, 
the 



2Ts 



-tt) sin(7r- 



Tg — to = to — Tyj which yields C(cos(^7r) — 
we obtain a maximum value ^ in [Ts,T]. 



■) = in [T^,Ts] giving the solution 
vr) sin(^7r)) (a minima). Similarly, 



t. 



Following the above procedure, it is easy to verify (from Eq. [T^ that at 
= T/2 the dispersion of GRF in the anteroposterior direction attains a minima 



ofA((^ 



''^'-■k) cosf^vr) 



sm 



Tin 



COSi 



^ v.v.^^ „ y ^...^ y^-7r)) and a maxima of i?((^^^7rj sml^^TTy --"V2rs 



vr)). 



in the mediolateral direction. Since, in gait, A > B > C, the sum of dispersions in 
three directions is minimal at to = |T. This confirms our results shown in Figs. HP 
- F and EH- Similarly, the dispersion of COM acceleration, dispersion of COM velocity 
and the COM mechanical energy consumption are all the minimal at to = |T. This 
phenomenon is independent of the physiological factors such as height, weight and gait 
parameters HI 



4. Conclusion 



In periodic motion, the foot's stance and swing substitute one another such that one 
foot always remains in stance. This allows the human body to be acted by the periodical 
GRF, which is related not only to the foot's movement style, but also to the substitution 
style of one foot with another. Following the biped movement style in normal gaits and 
the similar traits of GRF, we have been able to study the effect of tIFC upon the human 
COM dynamic characteristics based on the assumption that both feet's GRFs are the 
same in the same phase. Our results suggest that when tIFC falls in the middle of the 
other foot's stride cycle, the COM kinematics and dispersions of GRF acted on COM 
are the minimal, which has entailed the minimal average MEG of muscles. Our analysis 
suggests that it falls into the category of least-action principle and is consistent with 
the phenomenon that exists in normal gaits [T8] . 

Based upon the least-action principle, we have observed that tIFC has caused the 
GRF and the COM regularities to form a closed-curve on the transverse and frontal 
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plane, which might present a new method for gait evaluation. We advocate the use 
of this method to uncover and speed up the diagnosis simply by measuring the GRF 
for the patients with foot injuries or arthritis (for such patients, the model for these 
quantities is not symmetric [27]). Meanwhile, the patients need only walk a few steps in 
their normal gait, their COM dynamic characteristics will be acquired easily and more 
accurately. This is exactly what the clinic diagnosis is looking for. In addition to the 
human gait's adaptation to natural environment [25l [28] , the evolution of gait is the 
result of its observation of least-action principle. We believe that precise measurements 
of the variations of shear stress in natural gaits (bare-footed) shall profoundly enriched 
the content of least-action principle [29]. A further study of least-action principle will be 
significant to the domains such as sport rehabilitation, biometric identification [30| [31] 
and control of biped robot gaits [32l |33l [34] • 
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